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Abstract 

In this paper we present the exact solution of one dimensional Schrddinger equa- 
tion for Woods-Saxon plus Rosen-Morse plus symmetrical double well potential via 
Nikiforov-Uvarov mathematical method. The eigenvalues and eigenfunctions of this 
potential are obtained. The energy equations and the corresponding wave function 
for special cases of this potential are briefly discussed. The PT-symmetry and Her- 
miticity for this potential are also considered. 

Keywords: Schrddinger equation, Woods-Saxon, Rosen-Morse, symmetrical double 
well potential, Nikiforov-Uvarov 
PACs No. 03.65. Ge; 03.65.Pm 

1 Introduction 

It is well known that the exact solutions of the Schrddinger equation for some physical 
potential, plays an essential role in non-relativity quantum mechanics and are possible 
only for a few set of quantum systems. 

In recent time, a number of researchers have studied the solution of the Schrddinger 
equation. For example, these investigation have employed Woods-Saxon Potential [H El El S] , 
Harmonic oscillator potential [5], Hulthen potential [6j U\ E], Rosen-Morse potential [9J, HO] . 
Manning-Rosen potential [TT] etc. Several methods such as the Nikivorov-Uvarov method, 
Supper-Symmetric quantum mechanics, Factorization method, etc, have been used to solve 
differential equations arising from these considerations. 

In this present work, our focus is to solve the one dimensional Schrddinger equa- 
tion for Woods-Saxon plus Rosen-Morse plus symmetrical double well potential using the 
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Nikiforov-Uvarov mathematical method. The results include the PT-symmetric and the 
non PT-symmetric solution of this potential 

The paper is organized as follows, section 2; we briefly review Nikiforov-Uvarov (N — U) 
Mathematical method. In section 3, by using this method, we obtain the energy eigenvalue 
and eigenfunction for real and complex cases of this potential. We investigate the PT- 
symmetric and non PT-symmetric of this potential in sect. 4 and sect. 5. In section 6, we 
investigate the special cases of this potential and section 7 devoted to conclusion. 

2 The Nikiforov-Uvarov (N-U) Method 

The N-U method is based on solving a second order linear differential equation by reducing 
it to a generalized equation of hypergeometric type [121 dS|- By introducing an appropriate 
coordinate transformation z = z(r), this equation can be re-written in the following form. 

where a(z) and a(z) are polynomials, at most of second degree, and f(z) is a first degree 
polynomial. Now, if one takes the following factorization 

V(z) = y(z)<f>(z). (2) 

The equation (pQ) reduces to a hypergeometric type equation 

a{z)y"{z) + r(z)y'(z) + Xy(z) = 0, (3) 

where 

r(z) =t(z) +2tt(z). (4) 
The function n and the parameter A require for this method are define as follows 



a'(z) - t(z) 



and 



a'(z) — t(z)\ 2 



a 



z) + ka{z) (5) 



A = k + tt'(z) (6) 

In other to find the value of k, the expression under the square root must be a square 
of a polynomial. This gives the polynomial tc(z) which is dependent on the transformation 



function z(r). Also the parameter A defined in equation ([6]) takes the form 

n(n — 1) 



A = \ r 



-TIT 



-a 



The polynomial solutions y n (z) are given by the Rodrigue relation 

N d n 

y n ( z ) = — [a n {z) P {z)] 
y ' p(z) dz n 1 v >Hy n 

where N n is the normalization constant and p(z) is the weight function satisfying 

j- z W(z)p(z)] = r(z)p(z). 

Second part of the wave function cj){z) can be obtained from the relation 

d 

tt(z) = a(z) — [\n(j)(z)] 



(7) 



(8) 



(9) 



(10) 



3 Bound state solutions 



The one dimensional time independent Schrodinger Equation can be written as 
d 2 ^(x) 

■[E-V(x)]V{x) = Q , for h = 2p=l 



dx 2 



The generalized Woods-Saxon plus Rosen- Morse plus generalized symmetrical double well 
Potential can be expressed as 



V(x) 



l+ ge -2a* + (1+^-2^2 + ^3 SeC fl^OiX) 



\ +X4tan h q (ax) — V5 tan h 2 (ax) + Vq sec h 2 (ax) 
Using the following deformed hyperbolic function we have 



e — qe e + qe sin h q x 

sin h Q x = , cos h Q x = , tan h„x = : — 

2 2 cosh n x 



cos h n x ' 



cosech„x = 

sin h n x 



1 , cos h q x 

, COt h q X - 



sin h q x 



By using 



We obtained 



dx 2 



-lax 



d 2 d 
4a 2 z 2 -— + 4az 



dz 2 



dz 



(12) 



(13) 



(14) 



(15) 



3 



Putting Eqs. (115]) . f jl4j) . (fT2j) and appropriate parameter in Eq. ( TTBl) into Eq. ( TTT|) . we 

have 

(*) + n'^V ^ + 2^1 1 \2 {ft? - A - ? 2 (A + ft + 0} 

+ ^ {-A - 4 (ft + ft) - 2gft + 2^} + (ft - £)] (16) 
Where we have used the following notations 

^=4^ and ^ = h ' Z = 1 " 6 (17) 
Comparing Eq. (1161) with Eq. ([T]), we have 

a (z) = z 2 {ftg - ft - </ 2 (ft + ft + £)} + z {-ft - 4 (ft + ft) - 2gft + 2^} + (ft - £) 

f (z) = 1 — gz , er(z) = z — gz 2 (18) 
Inserting the values of these parameters with cr'(z) = 1 — 2qz into Eq. (J2]), we obtained 

tt(z) = ~Y ± \ K 25 - ^H} qz - 26] , for k = 2q (ft - ft) - 4 (ft + ft) - ft + V Hq6 

(19) 

Where we have used the following notations 



i? = ±l, 5= V / ^ft,^ = ^ + 16^±^ + 12ft + l (20) 

By choosing an appropriate parameter of k in ir(z) which satisfies the condition f(z) < 0, 
one gets 

„{z) = -^-^[{6-r l H}qz-28\ for k = 2q (ft - ft)-4 (ft + ft) -ft +7/^5 (21) 

From Eq. (jSJ), we obtained the parameter A as 

A = 2q (ft - ft) - 4 (ft + ft) - ft + V Hq6 - | - ~ {{25 - V H} qz} (22) 

Also we obtain the parameter A from Eq. ([7]) as 

A = A n = nq{(n + 1) +28-rjH)} (23) 
Comparing Eqs. (122|) and ( |23|) . we obtained the exact energy eigenvalues equation as 



n > 0,q > 1 



(24) 



Where we have substituted for /3i, /? 2 , 03, Pa, Pb, fle, 5 and H. From Eq. (jUJ) it 
can be shown that the weight function p(z) is 



= z 25 (l - qz) 



-T)H 



(25) 



Substituting Eq. fl25J) into the Rodrigue relation (jSJ), we obtained 



AT n P^-^)(l-2^) 



Where, we have used the properties of Jacobi Polynomial [151 
part of the wave function from Eq. Q as 

<f>(z) = z^l-qz)^ 1 '^ 



(26) 



We obtain the other 



(27) 



Finally we obtained the wave function \l/ n (z) 



z) y n lz) as 



' 1 7 



n 1 1" 



■2^) 



(28) 



Where iV n is the normalization factor to be determined from normalization condition 

f 1 \V n (z)\ 2 dz = l (29) 

JO 

This can further be written as 

N l f 1 z 2c (1 _ g ^(l+2d) [pM^)] 2 = j (30) 



Where c = 8 and d = — ^- . Using the different forms of Jacobi polynomial [T5| [T6] 

P M) ^ = 2 -n + C ) + (1 - ^(1 - zf 



Where 



r=0 

(rrz + c+ 1) 
n!r(n + c + d + 

_ F(n+1) 



r=0 



n \ r(n + c + d + r + 1) - 1 

r 



r!(n-r)! T(r+l)r(n-r+l) 



We obtained 



p(2c,2d) ( 



r + c + 1 



from Eqs. flS) and ([32]) as 



(31) 
(32) 



(33) 
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Putting Eq. (133]) into (130|) . we have 

N 2 n A n q (p, r) f 1 ^«+r-P+2«(i _ gz) 1+2d+p dz = 1 (34) 

V 

Before presenting the normalization constant N n , let us recall a relation of hypergeometric 
function, which is used to solve Eq. (1341) and present N n 



° af-^l - qz)- b ds = 2 F 1 (a, b;a + l;q)x g) (35) 

i r( c ) 

Provided Re(c) < Re(a) > 0, |arg(l — q)\ < n 

I z a -\l - q Z y b ds = -[ 2 Fx(a, b; a + 1; g)] (36) 

*/ CJ 

Using the above Eq. (I3"B"1) . we obtain the normalization factor N n as 



Where 



B nq (p,r) 
\| A nq (p, r) 



(37) 



(-l) n [r (n + 2c + l)] 2 r(n + 2d + 1) (-l) p+T g TO - p+r T(n + 2c + 2d + r + 1) 

n?tP» r ) - r ( n + 2c _ p + i) p (r + 2d + 1) T(2c + 2d + 1) ^ p\r\ (n - p)\ (n - r)W(p + 2d + 1) 

(38) 

(Bn^p.r))- 1 = — 1 ——— 2 F 1 (n + r-p + 2S+l 1 r]H-l-p;n + r-p + 2S+l;q) 

n + r — p + 2o + l 

(39) 



c=6, 5 = J£-fi A , f3 4 = —, d= — — , H = J—— + A — + 3— + 1 40 

v 4cr 2 V Q a Q a a 

4 PT -Symmetry and non-Hermiticity 

When a Hamiltonian commutes with PT, it is called PT-symmetric Hamiltonian, i.e. 
[H, PT] = 0, where the parity operator P and the time reversal operator T, satisfying the 
following relations. 

PpP' x = —P = TpT~\ PxP" 1 = -x, TAT' 1 = -%A (41) 

In this case, we set the potentials parameter in Eq. (Tl2|) as Vi, V2, V3, V4, V5, Vq, 
g 6 3 and a G flS (a — >■ za). Where S and nSf denotes the set of purely real numbers and 
purely complex numbers respectively. 
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This replacement and after some manipulations in Eq. ( Fl2|) leads to the following 
potential 



— Vi [1 + g(cos 2ax — isin2ax)] + V2 [cos 2ax — isin2ax] — 4(V3 + VJj) 
-V4 [(1 ~ g 2 )cos2ax + i(l + g 2 )sin2ax] + V 5 [(1 + q 2 ) cos2ax + i(l - g 2 )sin2ax - 2g] 



(1 + q 2 ) cos 2ax + (1 — q 2 ) isin2ax + 2q 



(42) 



Using Eqs. (iTTjl and ( 142]) and making the corresponding parameter replacements in Eqs. 
24|) and ( |28|) . we obtained energy eigenvalues equation and eigenfunction respectively as 



a 



Vx-lqVA-qVs _ V 2 
qd 1 q 2 a z 



2n + 1 - v Jl - 4^-4^1 -3^ 

'V q z a z qa* a z 



' 2n + 1 _J 1 _J^_ 4 gi+M-3 F5 



q 2 a 2 



qa" 



a' 



■Vi 



n > 0,g > 1 



(43) 



^ n (z) = CUV^(1 - ^) 



5) p^^-V 1 -*- 4 ^- 3 ^) (1 . 



-2^) 



(44) 

Using the above Eq. f|44|) . we can consequently obtain the normalization constant C n as 



Bng(p,r) 

\ Anq(P,r) 



Where 
A nq (p,r) 



(-l) n [T (n + 2c + l)] 2 r(n + 2d + 1) 



p=o 



(45) 



(-l) p+r q n ~P +r T(n + 2c+2d + r+l) 



r (n + 2c - p + 1) r (r + 2d + 1) T(2c + 2d + 1) ^ p!r! (n - p)\ (n - r)T(p + 2d + 1) 

(46) 



Bnq(P,r)) 1 



1 



n + r — p + 25 + 1 



2-Fi (n + r — p + 25 + 1, r\H — 1 — p; n + r — p + 25 + 1; q) 

(47) 



4a 2 



q 2 a 2 



Q- 



5 Non PT -symmetry and non Hermiticity 



We obtain another form of the potential by assuming some of the parameters as pure 
imaginaryi.e. q — > iq, a — > ia, V\ -> iVi, V3 — > iV$, V§ — > Wq. Under this replacement 
and after some manipulations, the potential takes the form 
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V{x) 



V\ [qcosax + i (1 + qcos2ax)] + V^cos 2ax — isin2ax] — i{V^ + Vq) 
-Va [(1 + q 2 ) cos lax + i(l - q 2 ) sin 2ax] + V 5 [(l + g 2 ) cos 2ax + z (l - q 2 ) sin 2ax - 2ig] 



(1 — q 2 ) cos 2ax + i(l + q 2 ) sin 2<xr + 2ig 



(48) 



Using Eqs. (II ip and (I48p and making the corresponding parameter replacements in Eqs. 
and ( 12"8"|) . we obtained energy eigenvalues equation and eigenfunction respectively as 



F - a 



qol 2 



q 2 a 2 



2n + l-^/l + ^ 



' 2n + 1 _J 1 + i_ 4 2i±M-35 



q 2 a 2 



qa" 



-F 4 



n > 0,g > 1 



(49) 



# n (z) = D n zVZ M(l - gz) 2 V V « 2q2 9 q2 a J P n K V 7 (i- 

(50) 

Using the above Eq. (j50p . we can consequently obtain the normalization constant D n 



■2qz) 



as 



Where 



Anq(P,r) 



D r 



Bnq(P,r) 

\ A m (p, r) 



(51) 



(-l) n [r(n + 2c+ l)] 2 r(n + 2d+ 1) (-l) p+r q n -P +r T(n + 2c + 2d + r + 1) 



r (n + 2c - p + 1) r (r + 2d + 1) T(2c + 2d + 



1) 4^ v\r\ 



p\r\ (n - p) \ (n - r)\T(p + 2d + 1) 

(52) 



^(p.r))" 1 



1 



n + r — p + 26 + 1 
V 4 



2F1 (n + r — p + 28 + 1, — 1 — p; n + r — p + 25 + 1; q) 

(53) 



4a 2 



q 2 a 2 



a- 



6 Discussion 



In this section, in the framework of the Schrodinger equation the energy eigenvalues and 
the corresponding wave functions for the Woods-Saxon potential, Rosen-Morse potential 
and symmetrical double well potential are obtained by setting appropriate parameters in 
Eq (I12p to zero. The PT-symmetric and the non PT-symmetric solution of this potential 
are also considered 
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6.1 Woods-Saxon Potential 



On putting V3 = V4 = V5 = Vq = 0, the potential (JT2]) turns to the Woods-Saxon 

pa ei si s 



-2ctx 



+ 



.1/ e -4ax 



1 + qe~ 2ax 1 + qe~ 2ax 
The energy eigenvalues equation and wave function are respectively obtained as 

2 



(54) 



E„ 



a 

T 



2n + 1 



q Q 2 q 2 ct 2 



V 2 



1 



2n + 1 - 



g 2 a; 2 



v 2 



in 



v 2 



n > 0, q > 1, 77 = 1 
(55) 
(56) 

This is consistent with [TB] 

Using the above Eq. (l56p . the normalization constant P n is obtained as 



P„ 



B nq (p,r) 

\] Anq (p, r) 



(57) 



Where 



(-l) n [r (n + 2c + l)] 2 r(n + 2d + 1) _ ^ (-l) p+ V~ p+r T(n + 2c + 2d + r + 1) 



r (n + 2c - p + 1) T (r + 2d + 1) T(2c + 2d + 1) ^ p\r\ (n - p)\ (n - r)\T(p + 2d+l) 

(58) ' 



[B nq (p,r)) 1 



n + r — p + 28+1 



iF\ (n + r — p + 28 + 1, t]H — 1 — p; n + r — p + 28 + 1; q) 



(59) 



c = 8 — d 



T]H 



H 



Vo 



2 ' " ^q 2 a 2 + l 

Now, let consider the case in which a is a pure imaginary parameter, i.e. a — > ia. Under 
this replacement and after some manipulations, the potential takes the form 

-Vi [1 + q (cos 2ax — isin2ax)} + V2 [cos 2ax — isin2ax] 



V (x) - 

(1 + q 2 ) cos2ax + (1 — q 2 ) isin2ax + 2q 
We obtain the real positive energy eigenvalues as 



(60) 



rp _ 



Vi . V2 



- 271+1 



V2 

q 2 a 2 



r] = l, n>0, q>l (61) 
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The above Eq. (IBTj) is consistent with [16] and [T7] for h = 2fi = 1 and a = By 
assuming some of the parameters as pure imaginary, i.e. q — > zg, a — > iot and Vi — > iV±, 
then we have potential of the form 

Vi [qcosax + + qcos2ax)] + V^cos 2ax — isin2ax] 



V (x) 

(1 — q 2 ) cos 2ax + i (1 + g 2 ) sin 2ax + 2ig 
We obtained the corresponding energy eigenvalues as 



E n — 

4 



Vi v 2 



2n + l 



v 2 



2n + 1 



' V 2 
q 2 a 2 



6.2 Rosen-Morse Potential 



(62) 



, n > 0, q > 1 (63) 



If we substitute for Vi = V2 = V5 = Vq — 0, then the potential in Eq. ( 112]) becomes 
Rosen-Morse 

V (x) = —V^sech 2 , (ax) — Vjtanhg (ax) (64) 

The energy eigenvalues equation and corresponding wave function are respectively obtained 

as 



E„ = — 



a 



2qV 4 
qo? 



2n+l + ,/l + 4^ 



2n+l + Wl + 4- 



-—\2n + l + 1/1 + 4 



ga 2 / a 2 



__1 

ga 5 



Ji 

ga 2 

-2 



- Vi 



n > 0, q > 1 (65) 



2 9 z), 77 = — 1 (66) 



Which is consistent with [TU] Using the above Eq. ( I66p . we obtained the normalization 



constant R n as 



Rr 



\ Kq (P, r) 



(67) 



Where 



(-l) n [r (n + 2c + l)] 2 r(n + 2d + 1) (-l) p+r g n " p+ T(n + 2c + 2d + r + 1) 

r (n + 2c - p + 1) T (r + 2d + 1) T(2c + 2d + 1) ^ p!r! (n - p)! (n - r)\T(p + 2d + 1) 



(■Bng(p,7-)) 1 



1 



n + r — p + 25 + 1 



(6* 

jFi (n + r — p + 25 + 1, r^i? — 1 — p; n + r — p + 25 + 1; q) 

(69) 



2 ' 



qa' 
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We consider a case in which a is a pure imaginary parameter, i.e. a — > lot. Then Eq. ([6 
takes the form 

-4V3 - I/4 [(1 - q 2 ) cos 2ax + % (1 + g 2 ) sin 2ax] 



V (x) 

(1 + q 2 ) cos2ax + (1 — q 2 ) isinlax + 2g 
The positive energy eigenvalues is given as 
■2 



(70) 



E n = — 



2\ 4 

«2 



2n + l- A /l-4^ 



2n+ 1 - 1 1-4- 



ii 



-Vi, n>0,q>l (71) 



Setting some of the parameters as pure imaginary, i.e. q — > iq ,a — > ia, V 3 — > iV 3 and 
V4 — > iV4, Then Eq. (|64|) takes the form 

— iV 3 — V4 [(1 + q 2 ) cos 2ax + z(l — q 2 ) sin 2ax] 



y (x) 

(1 — q 2 ) cos 2ax + i (1 + g 2 ) sin 2ax + 2iq 
The energy eigenvalues equation is also found to be (1TTT) 

6.3 Symmetrical Double Well Potential 

Choosing V x = V 2 = V 3 = V A = in Eq. (JT2} we have 

V (x) = V5 tanh 2 , (ax) — Vssech 2 {ax) 



(72) 



(73) 



Again, energy eigenvalues and the corresponding wave function are respectively obtained 

as 

,2 



En 



of 

"T 



14 

, v 2 



2^ + 1-^4^ + ^ + 1 



2n+l-W4-^ + ^ + l 



-/3 4 ,-r,,/4^% + ^+l 



, n>0, q>l 
(74) 

1 -2qz),rj = ±1 
(75) 



Using the above Eq. (1751) . we obtained the normalization constant W n as 
Where 



Bnq(p,r) 
\ A nq (P, r) 



(76) 



(-l) n [r (n + 2c + l)] 2 r(n + 2rf + 1) (-l) p+ y-P+T(n + 2c + 2rf + r + 1) 

r (n + 2c - p + 1) T (r + 2d + 1) T(2c + 2d + 1) ^ p!r! (n - p)! (n - r)\T(p + 2d+l) 

(77) 
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(B nq (p, r)) 1 = — ■ 1 2^1 (ri + r - p + 28 + l,rjH - 1 - p; n + r - p + 25 + 1; q) 

n + r — p + 2o + l 

, (78) 

c ft a „ F. v 6 3K 

c = 5 = J£, d= — — , H=J4 — 2- + ^ + l 
v 2 V Q a a 

Let us consider the case where at least one of the potential parameters is complex. If a is a 
pure imaginary parameter, i.e. a — > ia, symmetrical double well potential can be written 

as 

T/ . -4V 6 + V 5 [(l + q 2 )cos2ax + i(l-q 2 )sm2ax-2q] 

V (x) = ^r- ^- (79) 

(1 + q 2 ) cos2aa: + (1 — q 1 ) isin2ax + 2q 



By using Eqs. (QJJ) and (\78\i and making the corresponding parameter replacements in Eq. 
021]) again, we obtained the real positive energy eigenvalues equation 

.2 



TP - a ' 



qv 5 / , "1 2 



f y y - 2n+l- Jl- 4 A_ 3 ^ 



, n > 0,q > 1 
(80) 



When a — > ia, V G — > iV 6 and g,6 3 then V{x) takes the form 

-iV 6 + V 5 [(1 + q 2 ) cos 2ax + i (1 - q 2 ) sin 2ax - 2iq] 
(1 - q 2 ) cos 2ax + % (1 + g 2 ) sin 2ax + 2zg 



The positive energy eigenvalues is also found to be (I8"U|) 

7 Conclusion 

In this paper we solve the Schrodinger equation for Woods-Saxon plus Rosen-Morse plus 
Symmetrical double well potential via Nikiforov-Uvarov mathematical method. The results 
include energy spectrum of Wood-Saxon potential, Rosen-Morse potential and Symmet- 
rical double well potential which are obtained by setting some parameters to zero. The 
PT-symmetric and non-PT-symmetric of Wood-Saxon potential, Rosen-Morse potential 
symmetrical double well potential have also been discussed. 
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